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I. INTRODUCTION 

Measurements of weak boson pair production at hadron colliders are vital for testing 
the standard model and probing beyond it. These processes are particularly important 
because they can be used to probe the electroweak symmetry breaking mechanism and test 
the triple weak boson couplings [|T|. The process pp — > Z7 is of interest as a test of the 
standard model. In addition, this process is sensitive to contact interactions which appear 
in composite models of gauge bosons [Q], and it can also be used to search for evidence of 
anomalous ZZ'j and Z77 couplings [^]. In order to perform such tests it is imperative to 
have precise calculations of hadronic Z7 production. Futhermore, the calculations should 
include the leptonic decay of the Z boson since the Z boson is identified by its leptonic 
decay products. 

At leading order in QCD [(9(a°)], hadronic Zj production proceeds via the quark- 
antiquark annihilation subprocess qq — > Z'-f. This process was first calculated in Ref. [Q. 
Zj events can also be produced at leading order by the photon bremsstrahlung process 
which proceeds via subprocesses such as qg -^ Zq followed by photon bremsstrahlung from 
the final state quark. This process becomes significant only at supercollider center of mass 
energies. The production of Z'-f in association with one or two jets was calculated in Refs. 
and 0], respectively. The gluon fusion subprocess gg -^ Z7, which proceeds via a quark box 
loop, was calculated in Refs. H,^. Although this process is of order a^, it can be important 
at supercollider energies due to the large gluon luminosity; the gluon fusion cross section is 
15-30% as large as the qq — > Z7 Born cross section at supercollider energies |P|. The QCD 
radiative corrections to hadronic Z'-f production were calculated in Ref. 0] for the case of 
a real Z boson in the final state. 

In this paper, the calculation in Ref. [|10| is extended to include the leptonic decay of the 
Z boson. Since the Z boson is identified via its leptonic decay products, the inclusion of the 
leptonic decay in the calculation makes it much more useful for comparing with experimental 
data. For example, cuts can now be applied to the final state leptons, thus allowing one to 



more closely mimic the actual experimental conditions. 

To date, the CDF and DO experiments at the Fermilab Tevatron collider have collected 
12 [|ri|] and 6 |T^ Z7 events, respectively. Both experiments are expected to collect data 
corresponding to an integrated luminosity of approximately 100 pb~^ by the end of the 
current collider run (Run 1-B). This represents a five fold increase in statistics over the 
presently available data sample. Larger cross sections and luminosities are expected at the 
LHC. Thus it should soon be possible to use hadronic Z7 production to probe the standard 
model. 

The remainder of this paper is organized as follows. The formalism used in the calculation 
is briefly reviewed in Sec. II, numerical results for the Tevatron and LHC center of mass 
energies are given in Sec. Ill, and summary remarks are given in Sec. IV. 

II. FORMALISM 

The formalism used in the calculation is reviewed in this section. The discussion will 
be brief since the formalism has already been described in the literature. The calculation 



is done by using a Monte Carlo method for next-to-leading-order (NLO) calculations |T^ 



in combination with helicity amplitude methods |T^. Details on the NLO calculation of 



hadronic Z7 production for a real Z boson in the final state can be found in Ref. ||T^ and 
helicity amplitude calculations of Z7 production are described in Ref. 0. 

The calculation is done using the narrow width approximation for the decaying Z bo- 
son. This simplifies the calculation greatly for two reasons. First of all, it is possible to 
ignore Feynman diagrams in which the photon is radiated off one of the final state lepton 
lines without violating electromagnetic gauge invariance. [Radiative Z decay events can 
be suppressed by a suitable choice of cuts; see Sec. IIIB.] Secondly, in the narrow width 



approximation it is particularly easy to extend the NLO calculation of Ref. |T0[ to include 
the leptonic decay of the Z boson. 



The QCD radiative corrections to hadronic Z7 production were calculated in Ref. fTO 



for the case of a real Z boson in the final state. The present calculation extends the results 
of Ref. [|T0| to include the leptonic decay of the Z boson. The procedure for incorporating 
the leptonic decay Z —* i't^ into the NLO calculation of hadronic Z7 production is the 
same as that used in Ref. [|T5| for incorporating the leptonic decay W ^ iu into a NLO 
calculation of hadronic W'j production. A detailed discussion of the procedure is given in 
Ref. |15|. Basically, except for the virtual contribution, all the NLO contributions for real 



Z'y production have the form 

a 



da^^^iqq -^ Zi) = rfa^°™(gg ^ Z7) 



l + C,^(... 



(1) 



where 0"^°'''^ is the lowest order Born cross section, Cp = 4/3 is the quark-gluon vertex color 
factor, and a^ is the strong running coupling. Thus the leptonic decays can be incorporated 
by simply making the replacement 

da^°'''iqq ^ Z7) > da^^'^'iqq -^ Z-f -^ rt-f) (2) 



in the formulas of Ref. [^ for NLO real Z7 production. 

The simple replacement described in the previous paragraph does not hold for the virtual 
correction. For this contribution we use the virtual correction for a real on-shell Z boson 
which we subsequently decay ignoring spin correlations. Since the virtual corrections are 
small [about 1% as large as the Born cross section] and the effects of spin correlations are 
small, the overall result of ignoring spin correlations in the virtual correction is negligible 
compared to the 20% - 30% uncertainty from the parton distribution functions and the 
choice of the scale Q^. [Note that spin correlations are included everywhere in the calculation 
except in the virtual contribution.] 

III. NUMERICAL RESULTS 

Numerical results for NLO Z7 production at the Tevatron [pp collisions at -^i = 1.8 TeV] 
and the LHC [pp collisions at -y/i = 14 TeV] are presented in this section. This section begins 
with a brief description of the input parameters and acceptance cuts. 



A. Input Parameters 

The numerical results presented in this section were obtained using the two-loop expres- 
sion for ttg. The QCD scale Aqcd is specified for four fiavors of quarks by the choice of the 
parton distribution functions and is adjusted whenever a heavy quark threshold is crossed so 
that as is a continuous function of Q^. The heavy quark masses were taken to be m^ = 5 GeV 
and mt = 174 GeV [|T^]. The standard model parameters used in the numerical simulations 
are Mz = 91.173 GeV, Mw = 80.22 GeV, a{Mz) = 1/128, and sin^ 9„ = 1 - {M^/Mzf. 
These values are consistent with recent measurements at the CERN e^e~ collider LEP [ |17[| , 



the Stanford Linear Collider [|T^], the CERN SppS collider [|I9[, and the Fermilab Tevatron 
collider |^. The soft and collinear cutoff parameters used in the NLO Monte Carlo formal- 
ism (see Ref. [0) are fixed to 6s = 10~^ and 6c = 10^^. The parton subprocesses have been 



summed over u, d, s, and c quarks. The leptonic branching fractions and the total width of 
the Z boson are B{Z -^ e^e+) = B{Z — > /i^/U+) = 0.034 and Tz = 2.487 GeV, respectively. 
A single scale Q^ = M| , where Mzj is the invariant mass of the Z'y pair, has been used 
for the renormalization scale /i^ and the factorization scale M^ (see Ref. [0] for a defini- 



tion of these scales). The numerical results were obtained using the Martin- Roberts-Stirling 
(MRS) |21[] set A parton distribution functions with A4 = 230 MeV. These distribution 
functions have been fit to next-to-leading order in the MS (Modified Minimal Subtraction) 



scheme |22|, which is the scheme used for the present calculation. 



B. Cuts 

The cuts imposed in the numerical simulations are motivated by the finite acceptance 
and resolution of the detector. The finite acceptance of the detector is simulated by cuts on 
the four- vectors of the final state particles. These cuts include restrictions on the transverse 
momentum pj, and rapidity y of the photon and leptons. 

Since detectors are generally best equipped to measure photons which are isolated, i.e., 



not accompanied by a large amount of nearby energy p3[, a photon isolation cut will be 



applied in the numerical simulation. First, the leptons and photon are required to be 
separated in the pseudorapidity-azimuthal-angle plane 

21 1/2 



Ai?(£,7) 



A^j' + fAr^,^'' 



(3) 



A separation cut of Ai?(£, 7) > 0.7 will be imposed. In addition, a photon isolation cut 
typically requires the sum of the hadronic energy i?had in a cone of size Rq about the 
direction of the photon to be less than a fraction eh of the photon energy E^, i.e., 

J2 ^had < eh ^7 , (4) 

Ai?<_Ro 

with AR = [(A0)2 + (Ar7)2]V^- A photon isolation cut with eh = 0.15 Hll and Rq = 0.7 



will be applied in the numerical results presented in this section. This cut suppresses the 
portion of the cross section that comes from the photon bremsstrahlung process, which is 
beneficial since this process is not well understood theoretically. 

Before summarizing the cuts, it is useful to discuss cuts which will suppress radiative Z 
decays. Since photon radiation from the final state lepton lines is ignored in the calculation, it 
is necessary to impose cuts which will efficiently suppress contributions from these diagrams. 
In radiative Z decays the lepton-photon separation in AR space peaks sharply at small values 
due to the collinear singularity associated with the diagrams in which the photon is radiated 
from a final state lepton line. A separation cut on AR{i, 7), which has already been adopted, 
will therefore suppress radiative Z decay events. Contributions from radiative Z decay can 
be further reduced by a cut on the i~i~^'~f invariant mass. In radiative Z decays the i~i~^'y 
system has an invariant mass M{i'i'^'~f) close to M^, whereas for Z7 production M{i~i~^'y) 
is always larger than M^ if finite Z-width effects are ignored. Thus i'i'^'y events originating 
from radiative Z decays can be suppressed if M{i~i~^'y) is required to be slightly above M^. 
Note that this requirement is automatically satisfied in the present calculation since putting 
the Z boson on the mass shell and demanding the photon to have ^7^(7) > p™™ requires that 

M(r£+7) > p??'° + [m| + ip^'"'f] '^' . (5) 
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Thus no additional cuts are needed in the numerical simulation to suppress the radiative Z 
decays. 

The complete set of cuts can now be summarized as follows. 

Tevatron LHC 



pr^ici) > 10 GeV 
p^ii) > 20 GeV 

\yii)\ < 1-0 

\y{e)\ < 2.5 
AR{£,j) >0.7 



Pr(7) > 50 GeV 

p^ii) > 25 GeV 

|y(7)|<2.5 

\y{i)\ < 3.0 
AR{i,-f) >0.7 



E 



Afi<0.7 



Eh <0.15E^ 



C. Cross Sections 



The total LO and NLO cross sections for the process pp ^ Z'-f + X —>■ £~i~^'~f + X are 
given in Table 1 for center of mass energies corresponding to the present Tevatron [^/s = 
1.8 TeV], an upgraded Tevatron [y/s = 3.5 TeV], and the proposed LHC [y/s = 14 TeV]. 
The cross sections in Table 1 have been summed over £ = e, /i and include the cuts listed 
in Sec. IIIB. Note that the LHC cross sections are smaller than the Tevatron cross sections 
due to the more restrictive pj^ cuts used for the LHC case [in particular the ^^(7) cut]. 

The differential cross sections presented here have also been summed over i = e, fi and 
include the cuts listed in Sec. IIIB. The figures are arranged in two parts, with parts a) and 
b) being the results for the Tevatron [^/s = 1.8 TeV] and LHC [^/s = 14 TeV] energies, 
respectively. NLO and LO results are shown with solid and dashed lines, respectively. 

Figure |1| shows the distribution of the invariant mass of the i^i'^j system M{i~i~^'j). 
The NLO corrections are fairly uniform in the invariant mass, so that the NLO distribution 
can be approximately described by simply scaling up the LO distribution by a factor of 1.2 
(1.3) at the Tevatron (LHC). The NLO corrections are larger at the LHC due to the larger 
gluon density at the LHC energy. 



The rapidity distributions [in the laboratory frame] for the charged leptons ?/(£) and the Z 
boson yiZ^ are shown in Figs. and |^, respectively. [Both leptons have been histogrammed 
in Fig. ^] The NLO corrections are largest in the central rapidity region where the cross 
section is also the largest, however, the NLO cross section is uniformly enhanced over the 
LO cross section in the entire rapidity range. This fact is illustrated in Fig. ^ where the ratio 
of the NLO to LO cross section is plotted as a function of yijl)- The ratio has a constant 
value of 1.2 (1.3) at the Tevatron (LHC) center of mass energy. 

The angular distributions of the leptonic decay products contain information on the 
helicity of the Z boson and are simplest in the rest frame of the Z boson. Figure |^ shows 
the polar angle distribution of the t~ in the Z boson rest frame, measured with respect to the 
Z boson direction in the Z7 rest frame, i. e., cos 6'_ = p^- -p^ where p^- is the unit-normalized 
three-momentum of the t" in the Z boson rest frame and p^ is the unit-normalized three- 
momentum of the Z boson in the Z7 rest frame. If no cuts were imposed on the leptons, 
this distribution would exhibit a 1 + cos^ 6*. dependance, which is indicative of transversely 
polarized Z bosons. Unfortunately, the shapes of the distributions in Fig. ^ indicate that 
there is little or no polarization information left in the distributions. In the presence of cuts, 
the angular distribution of leptons produced via the decay of a gauge boson is, in general. 



dominated by kinematic effects rather than polarization effects [^. The sharp drops in the 
distributions near cos6'_ = ±1 are due to the kinematic cuts. The NLO corrections are once 
again uniform over the range of cos ^_ . 

In the distributions presented so far, the NLO corrections simply enhanced the LO 
distributions and produced little or no change in the shapes of the distributions. This is not 
the case for transverse momentum distributions; instead, the NLO corrections increase with 
the transverse momentum. The transverse momentum distributions of the charged leptons 
PtW, the Z boson pj.(Z), and the photon Pt{i), are shown in Figs. ^ 0, and |, respectively. 
[Both leptons have been histogrammed in Fig. ^.] In all three figures, the NLO corrections 
increase as the transverse momentum increases. The large corrections at high pj, are due to 
the contributions from the real emission subprocesses qq -^ Z'-fg, qg -^ Z'yq, and gq -^ Z'yq 
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which enter at order a^. 

Comparing the curves in Fig. ^ for the two different energies, one sees that the curves at 
the LHC energy do not turn over at small values of pj^li). This is simply due to the more 
restrictive cuts used for the LHC case. 

The NLO and LO prp{Z) distributions in Fig. ^have different behaviors at small values 
of prp[Z). At LO the Z boson and photon are back-to-back in the plane transverse to the 
beam direction, thus the Pt{Z) distribution cuts off at the value of the minimum Pj'(7) cut 
[Pril) = 10 (50) GeV for the Tevatron (LHC)]. The NLO cross section on the other hand 
has a finite but rapidly decreasing distribution for values oi prp{Z) below the minimum Pt{i) 
cut. In this region of small prp{Z), the photon is recoiling against a jet. 

In addition to showing the NLO and LO differential cross sections versus Pxil), Fig. ^ 
also shows the 0-jet and 1-jet components of the NLO inclusive cross section. Here a jet is 
defined as a final state quark or gluon with 

Py(j) > 10 GeV and \y{j)\ < 2.5 (6) 

at the Tevatron, and 

Pr(j)>50GeV and |y(j)l < 3 (7) 

at the LHC. The sum of the 0-jet and 1-jet cross sections is equal to the inclusive NLO cross 
section. The decomposition shows that the 1-jet component is small at the Tevatron energy, 
and is small at low values of p^il) but becomes the dominant component at large values 
of Pxil) at the LHC energy. Thus most of the NLO corrections at high p^i'l) are due to 
events with hard jets. This is also true for the p^i^) and Pt{Z) distributions. 

The increase of the NLO corrections with the transverse momentum is clearly illustrated 
in Fig. ^ where the ratio of the NLO to LO cross section is plotted as a function of Pxil)- 
The ratio increases from 1.2 to 1.5 (1.3 to 2.0) over the range of pj.(7) shown for the Tevatron 
(LHC) center of mass energy. 



IV. SUMMARY 

The QCD radiative corrections to hadronic Z7 production have been calculated to order 
Us with leptonic decays of the Z boson included. The inclusion of the leptonic decays makes 
the calculation more realistic since it is the leptonic decay products that are observed in an 
experiment. Distributions of the final state decay products have been given for the Tevatron 
and LHC center of mass energies. The calculation includes typical acceptance cuts on the 
final state leptons and photon. 

The calculation is done by using the Monte Carlo method for NLO calculations in combi- 
nation with helicity amplitude methods. With the Monte Carlo method it is easy to impose 
experimentally motivated acceptance cuts on the final state particles, also, it is possible to 
calculate the order Ug QCD corrections for exclusive channels, e.g., pp —^ Z'-f + jet. The 
narrow width approximation has been used for the decaying Z boson. This simplifies the 
calculation greatly since it is possible to ignore contributions from radiative decay diagrams 
without violating electromagnetic gauge invariance. Futhermore, in the narrow width ap- 
proximation it is particularly easy to extend the NLO calculation of real Z7 production to 
include the leptonic decay of the Z boson. Spin correlations are included everywhere in the 
calculation except in the virtual contribution where they can be safely neglected. 

In general, the QCD radiative corrections are uniform in invariant mass and angular 
distributions; these distributions are scaled up in magnitude by the corrections, but undergo 
little change in shape. In contrast, the QCD radiative corrections increase with transverse 
momentum, and as a result, the p^ distributions are significantly enhanced at high prp. The 
large corrections at high pj, are due to contributions from the real emission processes which 
enter at order a^. The QCD radiative corrections increase with the center of mass energy 
due to the increasing gluon density in the proton. These behaviors of the NLO corrections 
are qualitatively the same as those found in the NLO corrections to W^'-f [1^, ZZ, W'^W'^, 
and W^Z production p6|. 
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TABLES 
TABLE L Total LO and NLO cross sections for the process pp -^ Z7 + X ^ i^i'^j + X, 
for center of mass energies corresponding to the present Tevatron [^/s = 1.8 TeV], an upgraded 
Tevatron [\/s = 3.5 TeV], and the proposed LHC [\/s = 14 TeV]. The cross sections have been 
summed over I = e, fi and the cuts listed in Sec. IIIB have been imposed. 

^ (TeV) pp a (fb) 

1.8 pp LO 360 

NLO 420 

3.5 pp LO 540 

NLO 620 

14. pp LO 230 

NLO 310 
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FIGURES 
FIG. 1. Invariant mass distribution of the i~i~^'j system 

for the process pp — > Zj + X — > i~£~^'j + X. Parts a) and b) are for the Tevatron and LHC 

center of mass energies, respectively. The NLO (sohd Une) and LO (dashed hne) cross sections are 

shown. The cross sections have been summed over i = e, n and the cuts listed in Sec. IIIB have 

been imposed. 

FIG. 2. Same as Fig. |l| but for the rapidity distribution of the leptons. 

FIG. 3. Same as Fig. || but for the rapidity distribution of the Z boson. 

FIG. 4. The ratio [da^^'^ /dy{i)]/[da^'^ /dy{i)] plotted versus y{i) for the process 
pp -^ Zj + X ^ i^i^j + X. Parts a) and b) are for the Tevatron and LHC center of mass 
energies, respectively. The cuts listed in Sec. IIIB have been imposed. 

FIG. 5. Same as Fig. Q but for the angular distribution of the negatively charged lepton. The 
angle 9- is measured in the Z boson rest frame with respect to the Z boson direction in the Zj 
rest frame. 

FIG. 6. Same as Fig. |l| but for the transverse momentum distribution of the charged leptons. 

FIG. 7. Same as Fig. || but for the transverse momentum distribution of the Z boson. 

FIG. 8. Same as Fig. |l| but for the transverse momentum distribution of the photon. In 
addition, the 0-jet (dotted line) and 1-jet (dot-dashed line) cross sections are also shown. 

FIG. 9. The ratio [da^^'^ /dpj,{'y)]/[da^'~' /dprpi^j)] plotted versus ^7^(7) for the process 
pp -^ Z'j + X ^ i~i'^j + X. Parts a) and b) are for the Tevatron and LHC center of mass 
energies, respectively. The cuts listed in Sec. IIIB have been imposed. 
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